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1. Summary of Major Research Findings 

The goal of our research over the past three years has been to contribute 
to the development of a definitive global climatology of the surface heat 
balance over the world ocean using a combination of satellite and conventional 
observed data. After careful development and testing of a method for esti- 
mating climatological sensible and latent heat fluxes (Appendix A) , we con- 
structed an atlas of the components of the surface heat balance over the 
global ice-free oceans (Appendix B) which may be regarded as updating the 
work of Budyko (1965) and his colleagues. During the course of our research, 
we reached a number of conclusions about the role that satellite data may 
play in improving the estimates of heat balance components in the near future. 
Our heat balance results and recommendations for the use of satellite data 
are briefly summari-ed below. 

1.1 Estimates of Climatological Latent and Sensible Heat Fluxes 

In the paper by Esbensen and Reynolds (1981) we have demonstrated that 
is is possible to estimate long-term monthly mean latent and sensible heat 
fluxes over the ocean to within the <20% relative accuracy of the bulk aero- 
dynamic formulas, by using observations of the monthly mean surface wind 
speed and the monthly mean sea-air temperature and humidity differences. 

Even more remarkable is that it appears to be possible to make a useful 
estimate of the fluxes on a month- to- month basis from monthly averaged sur- 
face data. Both of these findings have important implications for the use 
of satellite data (see Section 1.5). 
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1.2 An Atlas of the Sunfaae Heat Balance 

Using a climatology of surface marine meteorological observations pre- 
pared by the National Climatic Center and the cloudiness atlas of Berliand 
and Strokina (1980) we have produced an atlas of the components of the sur- 
face heat balance over the global ice-free oceans (Esbensen and Kushnir, 

1981) . The atlas includes maps of the monthly mean distributions of the 
calculated heat budget components and of the basic variables used in the 
calculation. All the fields have a uniform spatial resolution and are pre- 
sented on a 4® X 5° grid. Although the resolution of our data set is coarse 
we find good agreement between our atlas results and several recent regional 
studies. Our estimates have the advantages of being global in extent and 
achieving approximate annual heat balance over the global oceans. 

The atlas has already been useful in the development and analysis of a 
coarse resolution global ocean model at the Climatic Research Institute and 
we e.xpect that it will continue to be useful as the institute staff begins 
to explore coupled ocean- atmosphere models. The gridded data set is avail- 
able to other investigators in a very portable format on magnetic tape. 

1.5 Recommendations fan the Use or Sazellite Da^a 

We believe that it is now possible to clearly define the role of satel- 
lite data in improving climatological estimates of the components of the 
surface heat balance and in monitoring their interannual variability. We 


recommend the following: 


a) High priority should be given to monitoring the monthly mean of the 
magnitude of the surface wind stress using satellite data. We 
believe that our results show that this information can be used 
directly in bulk aerodynamic formulas to determine the monthly mean 
sensible and latent heat fluxes to well within a relative error of 
■v20% wherever those fluxes are significant. 

b) High priority should be given to developing and deploying a robust 
buoy-satellite instrument system which measures three quantities - 
surface air temperature and humidity, and the sea surface tempera- 
ture. This is particularly important over regions of the global 
oceans where ships-of-opportunity provide inadequate coverage and 
the associated space- time sampling errors can approach 50% (Weare 
and Strub, lOSll . 

cj High priority should be given to obtaining the monthly averaged 

surface radiation balance over the global oceans by fully utilizing 
existing satellite data. We believe that climatological estimates 
can be obtained by combining existing radiative flux data from 
satellites with surface meteorological data from ships and buoys. 

In this context we suggest that cloudiness data be viewed more as 
a by-product of the analysis rather than a primary input variable 
to empirical radiation formulas. (See Appendix B for e.xamples of 
this type of formula.) Of the three recommendations, we believe 
that this step would be the most likely to lead to a significant 
quantitative improvement in the estimates of the net downward heat 


flux into the oceans . 
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2. Summary of Research Activities 

The research summarized here was performed in collaboration with Dr. 

R. W. Reynolds^ and graduate students Messrs. Y. Kushnir and E. I. Tollerud. 

We have received good support from the staff of the computing facility at 
the OSU Climatic Research Institute (CRI) and excellent secretarial and 
clerical support through CRI and the Department of Atmospheric Sciences. 

In addition to the research activities at CRI, the results of our work 
have been presented at several scientific meetings including the Symposium 
on Empirical and Model-Assisted Diagnosis of Climate and Climate Change, 

1979, Tbilisi, U.S.S.R., the Third Conference on Ocean- Atmosphere Inter- 
action, 1980, Los Angeles, CA, and the Workshop on the Application of Existing 
Satellite Data to the Study of Ocean Surface Energetics, 1980, Madison, WT. 

At the Madison workshop, the Principal Investigator co-chaired the Working 
Group on Wind Stress and Heat Fluxes. 

5. Continuing Research 

The heat balance computations obtained under this grant are being ana- 
lyzed and extended as part of the Climatic Research Institute's ongoing study 
of the dynamics of seasonal climate change. By combining the values of evap- 
oration implicit in the latent flux estimates from the heat balance atlas 
shown in Appendix B with the precipitation data of Jaeger (1976), we have 
obtained climatological estimates of evaporation minus precipitation over 


^Present affiliation: Climatic Analysis Center, NMC, NWS, NOAA, Washington, 

D.C. 20255. 
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the global oceans for the twelve calendar months. These estimates have been 
very useful in the analysis of water vapor transport in the OSU Atmospheric 
General Circulation Model. We are also proposing to extend the surface heat 
balance computation to land and ice covered areas of the earth under a con- 
tinuing grant from the National Science Foundation since budgetary constraints 
within NASA have forced the termination of this grant. 


4. Publications of Grant NSG 5508 

5. K. Esbensen and R. W. Reynolds, 1981: Estimating monthly averaged air- 

sea transfers of heat and momentiim using the bulk aerodynamic method. 

J. Phys. Ooeanogz'.j 11, 457-465. (See Appendix A.) 

S. K. Esbensen and Y. Kushnir, 1981: The heat budget of the global ocean: 

an atlas based on estimates from surface marine observations. Report No. 
29 Climatic Research Institute, Oregon State University, Corvallis, 
Oregon, 27 pp. and 188 charts (in press). (See Appendix B.) 
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ABSTRACT 

Air*t«a transfers of sensible heat, latent heat and owmentum arc computed from 25 years of middle* 
latitude and subtropical ocean weather ship da^a in the North Atlantic and Nonh Pacific usins the bulk 
aerodynamic method. The results show thm monthly averaied wind speeds, temperatures and huinidi- 
ties can be used to estimate the monthly averaavd sensible and latent heat ftuxes from the bulk aero* 
dynamic equations to within a relative error of -10%. The estimates of monthly averaged wind stress 
under the assumption of neutral stability are shown to be within of the monthly averaged non* 
neutral values. 


1. Introduction 

Estimates of the transfers of sensible heat, latent 
heat and momentum between the atmosphere and 
the oceans may provide important insight into the 
dynamics of interannual and decadal climate fluc- 
tuations. Sensible and latent fluxei combine with 
radiative fluxes to determine the major part of the 
net air-sea exchange. Momentum transfers play an 
important role in driving the ocean general circula- 
tion and particularly in determining the vertical 
structure of the upper ocean. 

The first systematic evaluation of the sensible and 
latent heat transfers over the world ocean was per- 
formed by Budyko (1963). He computed long-term 
monthly mean values of the near-suiface wind speed 
and of the air-sea temperature and humidity differ- 
ences, and evaluated the fluxes by the well-known 
bulk aerodynamic formulas with constant transfer 
coefficients (see Section 2). 

Budyko' s method depends on the validity of the 
following two conditions: 

1) The covanance of the wind speed with either 
the air-sea temperature difference or the air-sea 
humidity difference must be less than the product of 
their respective mean values (except for cases in 
which the fluxes are small). This condition was found 
to be satisfied in several important climatic regions 
of the world oceans for limited time periods by 
Malkus (1962), Budyko and Gandin (19^), Kraus 
and Morrison (1966) and Fissel et al. (1977). 


' Present iflUixtioa; Oimaiic Analysis Center, NMC, NWS. 
NOAA. Washington, DC. 20233. 


2) The transfer coefficients for sensible heat and 
evaporation must not be strongly dependent on 
either wind speed or air-sea temperature and humid- 
ity difference.^ for typical conditions over the 
oceans. Budyko (1974) concludes that the stability 
of the surface boundary layer, as parameterized by 
the bulk Richardson number, is almost neutral ex- 
cept in weak wind conditions. As we will show in 
subsequent sections, computations using recent 
developments in surface boundary-layer theory 
support Budyko's conch sion. 

Tlie first systematic evaluation of air-sea. momen- 
tum transfers on a global basis was presented by 
Hellerman (1967). His computations were based on 
the bulk aerodynamic formulas with a drag coeffi- 
cient for neutr^ surface boundary-layer conditions 
and wind data consisting of monthly-averaged wind 
roses with both wind direction and speed categories. 
The use of only one mean wind speed in each direc- 
tion category of a monthly averaged wind rose leads 
to underestimates of the computed monthly averaged 
wind stress by 10-30% (see Hellerman, 1965). 

Recently there has been considerable interest in 
updating the classical estimates of Budyko (1963) 
and Hellerman (1967). An excellent example of this 
type of work is Bunker's (1976) estimates of the sur- 
f^e heat balance over the Atlantic Ocean. The 
approach used by Bunker differs from the classical 
approach since he computed the turbulent fluxes for 
each simultaneous measurement of surface meteoro- 
logical variables. The computed fluxes within a given 
month are then summed, and finally are divided by 
the total number of observations. While data proc- 
essing technology has made Bunker's approach 


00r-367(V81/O4O4J7.09$O6.25 
C 1961 American Meteorological Society 
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feasible, more than 30 million ship reports would 
b.sve to be processed at the present time to obtain 
global estimates of the turbulent heat and momen* 
turn transports. The classical method clearly has a 
number of practical advantages which will be dis> 
cussed further in our final section. 

The main purpose of this paper is to provide a 
basis for determining to what extent the classical 
method of estimating turbulent fluxes of sensible 
and latent heat from monthly averaged surface 
meteorological variables can be used for studies of 
interannual and decadal climatic variability. A 
unique feature of our study is the systematic ex- 
ploration of the effects of boundary-layer stability 
and of the physics of the air-sea interface on the 
estimates of the sensible and latent heat fluxes using 
recent developments in boundary-layer theory. We 
also will discuss the importance of subility effects 
for the estimates of air-sea momentum transfers. 
Compuutions will be based on data taken from 
ocean weather ships in ihe North Pacific and North 
Atlantic Oceans since these data are the best avail- 
able long-term records of surface marine meteoro- 
logical variables. In the final section, we will discuss 
the implications of cur results for studies of climatic 
variability. 

2. The bulk aerodynamic calculations 

The well-known bulk aerodynamic formulas for 
the sensible heat flux H and the latent heat flux LE 
at the air-sea interface are 

H - pc^CnVaiT, - T,), (1) 

LE * pLCcV,(q, - q«). (2) 

where the variables V',, T« and c« are the wind speed, 
temperature, and specific humidity of the air at a given 
height (:,) which is within the surface boundary 
layer. The variables T, and q, are the sea surface 
temperature and the saturation specific humidity of 
air with temperature T, and sea-level pressure p,.- In 
this study we have used z« « 10 m and approxi- 
mated the potential temperature difference (8. - 6,) 
by (7, - 7«). The variable p is the average air density 
in the layer from the surface to is the specific 
heat of air at constant pressure, and L is the latent 
heat of evaporation of water. The quantities Cm and 
Cf are the so^alled transfer coefficients for sensible 
heat and evaporation, respectively. 

To estimate the fluxes from Eqs. (1) and (2), the 
transfer coefficients must be determined. Although 
these coefficients may be determined in a com- 
pletely empincal manner using regression analysis, 
boundary-layer similarity theory provides a ra- 


’ The varubtes V., T„ T,, 4, amlv. arc all aatumed to be aver- 
aies over a time-scale at least aa lon( aa the lifeume of bound- 
ary-tayer convective elententa. 


tional, physically motivated way to minimize the 
dependence of and Cg on empiricism. In par- 
ticular, similarity theory allows ’•■s inclusion of the 
eflfects of the stability of the atmospheric surface 
boundary layer, the measure of stability being the 
Monin-Obukhov length L which is defined as 



where r is the magnitude of the surface wind stress 
(see Section 7), Tf it the virtual temperature, and 
g is the gravitational acceleration. Here the overbar 
indicates the Reynolds averaging of the product 
v’Tr', where w is the vertical air velocity. Similarity 
theory also allows the effects of the rou^ness of the 
sua surface to be parameterized by a roughness 
length 

We have chosen to evaluate the truisfer coeffi- 
cients using the method proposed by Liuer al. (1979) 
with the exception that we have neglected the dif- 
ference between the “observed" sea surface tem- 
perature and the skin temperature of the water at 
the air-sea interface. The effects of boundary layer 
stability are included using the Businger-Dyer model 
(see Businger et al., 1971, Dyer and Hicks, 1970). 
.Also, the Liu ei al. model determines the roughness 
length :« and the corresponding lengths for the pro- 
files of temperature and specific humidity, Zt and 
z„; the lengths are not assumed to be equal. There 
is no free convection limit for unstable conditions 
as considered by Deardorff (1972); however, calm 
conditions with large unsuble air-sea temperature 
and humidity differences are rare over the oceans. 

The dependence of Cm on wind speed and air-sea 
temperature differences computed by the Liu et al. 
method are shown in Fig. I for a relative humidity 
of 70%. For any given wind speed and temperature 
difference, C« and Cg differ oy less than 10% (see 
Liu et al.). However, at wind speeds < 10 m s*' 
both Cm and Cf are strongly dependent on both wind 
speed and air-sea temperature differences. At wind 
speeds > 10 m $*', the stability of the boundary 
layer has very little effect on the transfer coeffi- 
cients and the dependence on wind speed is weak. 

The transfer coefficients proposed by Liu et al. 
(1979) agree well with those of Kondo (1975), and 
give a reasonable fit to the currently available data 
from the surface boundary layers of the ocean and 
atmosphere. As we will show in Sections 4 and 5, 
the conclusions of this study do not depend crucially 
on the choice of transfer coefficients. Although this 
may seem paradoxical to the reader given the im- 
portance of the transfer coefficients in obtaining 
reliable estimates of the fluxes, this study is con- 
cerned with the methodology of estimating monthly 
averaged and long-term climatological fluxes from 
the bulk aerodynamic equations. We believe our 
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Trontltf Ce«ff)c<*ni for S«mibi« H«ot 



FiC. 1. Th« transfer cocfKcicni for sensible beet from Liu er oL 
(1979) as a function of wind speed at 10 m above the sea surface 
for various air-sea temperature differences {XT • T, - f,). 
The neutral surface boundary layer curve is sJso indicated. All 
curves assume a boundary-layer relative humidity of 70%. 

conclusions concerning the method for obtaining 
flux estimates are valid for any physically reason* 
able choice of transfer coeflicients. 


( ) - 1V"‘ i ( )e. (6) 

••I 

and N is the number of simultaneous measurements 
of Tt, T, and V,. It should be noted that p, c, 
and L are not included in the averaging operation 
since these variations either play a negligible role 
in the surface flux estimates or can be incorporated 
using linear combinations of fi and (e.g.. Brook. 

1978). 

The classical method (see Budyko, 1963) uses the 
expressions 

ft - pc,C„V,(t. - t.h (7) 

LE - pLCtV,^. - 4 -). ( 8 ) 

where and Cg are obuined by substituting the 
sample means off., T„ q, and V, into the a^rithm 

for Ch »nd Cf. 

If measurements of surface meteoroiogical vari- 
ables are always simultaneous and dense in space 
and time, and if the resources for data processing 
are adequate, the sampling method provides the best 
estimate of the monthly averaged flux. However, 
for non-simultaneous and sparse dau sets, the clas- 
sical method may have botSi practical and theoreti- 
cal advanuges as discussed in Section 8. 


3. Two estimation methods for monthly aver:.„.;d 
heat fluxes 

In this section we distinguish between the sam- 
pling method and tht classical method for computing 
monthly averaged transfers of sensible and latent 
hear at the sea surface from the bulk aerodynamic 
equations. !n the sampling method, fluxes are 
estimated from the sample means of H and LE in 
Eqs. (1) and (2) for a given calendar month and 
geographic area. These sample means are defined as 

N ■ pc,CmV,{T, - T,), (4) 

LE - pLCtV,(q, - qj, (5) 

where the averaging operation defined by the over- 
bar is 


4. Data 

To quantitatively test the differences between 
the two methods, we have analyzed three-hourly 
marine meteorological observations from nine ocean 
weather ships. The study was limited to the 25-year 
period from 1948 to 1972 and the ships were selected 
to give good geographical coverage. Table 1 shows 
the ship positions and periods for which observa- 
tions were available. The position of each ship also 
is shown in Fig. 2 with the ocean surface current 
regions from Sverdrup et al. (1942). As can be seen, 
the ships have sampled both eastern and western 
middle-latitude oceanic regions including thr strong 
western boundary currents of the Kuroshio and 
Gulf Stream. Six of the ships have nearly continu- 
ous records for the entire 25-year period. 



Ftc. 2. Fosiuons of ocean weather ships used in thu study supenmposed o<* :;>i(facc currems (ram Svcrdnip ti al. (19C). 
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Taele 1 . Oc«*« wMlber *hip tutkm 4«u ttv4M. 


Sution 

Latitude 

Loogitiidc 

Ptriod studied 

A 


33^ 


B 

j4W*N 

51*W 

l»4t-72 

D 

44^ 

4I*W 

1930-72 

H 

36*N 

Ttrw 

1990-54 

j 

32Mi*N 

2trw 

1931-71 


30*N 

140^ 

l94»-72 

p 

3074 

145^ 

1990-72 

V 

34*N 

164*E 

1996-71 

X 

3TN 

153^ 

194t-53 


We processed the dau to eliminate wy observa- 
tions that were obviously in error or in which the 
weather ship was not wnthin 1* of the position pven 
in Table 1. We also eliminated dau when simul- 
Uneous observations of wind speed, wind direction, 
sea level pressure, and air, sea and dew-point tem- 
perature were not available. 


Heat riuMS (W/m*l *S' 4u.,' 



5. l^ng-tenn cllmatolof kal estimates 


Fig. 3b. A* in Fig. 3« e*c«pi for July. 


In this section we see how closely the classical 
method approximates the sample mean of rtux esti- 
mates from temporally dense observations taxen 
within a given calendar month for a period of 10 to 
20 vears. The principal result is shown in Fig. 3 for 
January and July. The classical method agrees with 
the method of sample means to within a few percent 
in both summer and winter and in all climatological 
regimes of the subtropical and middle-latitude oewns 
which we have investigated. The only exceptions 


to this conclusion occur when the fluxes are small. 
However, when the fluxes are small, percentage 
comparisons are meaningless and Fig. 3 shows the 
agreement to be within a few Watts per sejuare meter. 
From result? not shown here, we reached the same 
conclusions using only nighttime obseiwations. 
Thus, diurnal effects, whether they arc real or instru- 
mental do not lead -o significant differtnees between 
the classical method estimates and the sampling 


Htot Flurts (W/m*) for Jonuory 



Fio, 3». CUniElotogical »«a»ibW aad Ituat h«»i flu» «***"^ 
tor jMuury for «ch »hip compuwb by the 
Ubteisu. EM) the clESMcai method (ordiiuie). Th* dished lines 
indifTte relEtive errors of 10%. See test for detaus. 


method estimates. 

The authors believe the above conclusions are not 
likely to be changed by any physically reasonable 
choice of transfer coefficients. To demonstrate this 
point, the fluxes shown in Fig. 3 were recomputed 
with the assumption that sensible heat, water vapor 
and momentum are exchanged between the occm 
and the atmosphere in the same manner. (Thus, in 
the Uu ef aJ. (1979) algorithm, :r - - ^1.^* 

shown in Fig. 4, the transfer coefficient .or sensible 
heat using this unrealistic assumption has an extreme 
dependence on wind speed. The curves for Cf are 
similar. The lieet fluxes for January with this as- 
sumption are shown in Fig. 3. Again, the discrepan- 
cies between the classical method and the sample 
mean are only a few percent for the lateiit heat 
fluxes and about or less for the sensible heat 
fluxes. This result implies that the difference be- 
tween Budyko’s (1%3) flux estimates and Bunker's 
(1976) estimates for the Atlantic Ocean are due 
primarily to daU availability, dau reliability, and 
the use of different transfer coefficients. The differ- 
ences are not due to the fact that Budyko chose 
the classical method and Bunker chose the sampling 
method. 

.A formal explanation for the agreement between 
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Tront»«r Co*ffici.nt for Stniitri* Htot H«ot Fluiu (WAn*) for Jonuory 



the two methods of estimating the fluxes is obtained 
from an inspection of Table 2. First, as discussed 
in previous studies with cons tant tra nsfer coeffi- 
cients, the contribution of the Va'lT' and Va'lq’ 
terms to the total flux is more than order of magni- 
tude less than the contribution of the *Jid 

terms except when the total fluxes arc small. 
Second, the dependence of the transfer coeffleients 
on the wind speed and air-sea temperature and 
humidity differences is small enough so that the re- 
maining Reynolds-averaged terms involving Cm' and 
Ce' also are more than an order of magnitude less 
than the total fluxes. (Although not shown here, 
results from compuutions of the fluxes with lo 
- zt - Ic^d to similar conclusions.) 


Fig. 5. Climatological sensible and latent heat fluxes as in Fig. 3a 
but computed with transfer coeffleients for which » Zr “ c,. 

Finally, we consider the importance of including 
boundary-layer stability effects on long-term 
climatological estimates of the sensible and latent 
heat fluxes. The heat fluxes with the neutral trans- 
fer coefficients from Fig. 2 were recomputed using 
the sampling method, and Fig. 6 shows a compari- 
son of these results with the stability-corrected 
sampling method results. The figure shows differ- 
ences of ~10%. This result supports Budyko's 
(1974) conclusion which was reached using a bulk 
Richardson number as the stability parameter for 
the boundary layer. 


Table 2. Terms in the Reynolds-averaged bulk aerodynamic equations for boundary-layer temperature and moisture fluxes 
at representative ship locations. The transfer coeffleients of Liu tt al. have been used. 


Temperature flux (10*’ m s*' K) 



^ CmV.AT - 


C„v,z.r + 

V.ChXT- 


+ ChV.XT' 

Ship D 

Jan 

66,12 

62.86 

2.30 

3.28 

-3.35 

-1.17 

Jul 

2.45 

1.88 

-1.03 

2.44 

0.01 

-0.85 

ShipN 

Jan 

12.00 

10.81 

-0.36 

2.40 

-0.28 

-0.37 

Jul 

8.37 

6.30 

0.45 

Moisture flux (10~* m $*' 

1.32 

gg-i) 

0.03 

-0.16 


~ 


+ CfV,',3q' + 



-I- C.-V'.-Oq' 

Ship D 

Jan 

68.0 

70.7 

-0.1 

1.9 

-3,7 

-0.8 

Jul 

16.4 

16.5 

-1.9 

2.6 

-0.0 

-0.8 

Ship N 

Jan 

37.6 

39.1 

-1.4 

1.2 

-1.1 

-0.2 

Jui 

35.3 

35.2 

-0.0 

0.2 

-O.I 

-0.0 
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Haot FIu»m (W/m*) tor Jonuory 



Fic. 6. Clinwiological sensible and latent heat fluxes as in 
Fig. 3a except that valttes along the ordinate are computed by the 
sampling method with neutral transfer coefficients. 


6. Monthly averaged estimates 

The results from the previous section suggest that 
the classical method also might be successful in 
monitoring the interannual variability of the sensible 
and latent heat fluxes. To examine this oossibility we 
calculated the heat fluxes for each January and July 


Heat Fluxes for January for Ship 0 



Fig. 7s. Monthly averaged sensible and latent heat fluxes 
computed by the sampling and classical methods for each Janu- 
ary at Ship 0. 


Heot Fluxes for Jonuory for Ship N 



at all the ocean weather ships using both the classi- 
cal method and the sampling method. Since the re- 
sults of the comparison are almost independent of 
the choice of month or ship, we have chosen to 
present only the January results at ships D and N 
representing significantly different heat flux regimes. 

Fig. 7a shows the comparison for Ship D. Although 
the average heat fluxes are large (see Fig. 3), (‘'e 
interanntial variability of the monthly averaged 
fluxes is also large. Despite the observed variability, 
the classical method still gives estimates of the heat 
fluxes which are within 10% of the values obtained 
by the sampling method. These results are consistent 
with the conclusions of Fissel et al. (1977) for 
Ship P; thev performed cross-spectral analysis on a 
series of 2-year segments of data at Ship P with con- 
stant transfer coefficients and found that disturb- 
ances with perioos of less than a month contributed 
onlv 10% of the biennially averaged sensible heat 
flux and only 5% of the biennially avcv.jged latent 
heat flux. 

The results in Fig. 7b for Ship N show that the 
interannual variability of the heat fluxes in subtropi- 
cal regions also is well resolved by the classical 
method. The classical estimate of the latent heat 
flux is systematically higher than the sampling esti- 
mate. but still within 10%. 

The comparisons in this section were based on the 
transfer coefficients shown in Fig. 1. However, we 
also computed the monthly averaged fluxes with the 
unrealistic transfer coefficients shown in Fig. 4. 
Since these coefficients have a stionger dependence 
on wind speed, the sampling and classical methods 
arc not in such close agreement. However, the average 













w^im^wfjm jwij 
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disagreement for the 25-year period was still < 10%, 
and never exceeded 20% in any individual month. 

Finally, we wish to point out that the transfer coef- 
ficients themselves are subject to uncertainties of 
-20-30% (e.g.. Pond et at., 1974). In view of this 
uncertainty, the classical method is a practical and 
reliable alternative to the method of sample means. 
The uncertainty underscores the need for independ- 
ent methods for obtaining the heat fluxes such as the 
budget approach of Oort and Yonder Haar (1976). 

7. CUmatoiogical ssd monthly averaged wind-stress 
estimates 

As discussed by HeUerman (1965) a wind rose 
with both direction and speed categories is neces- 
sary to accurately determine the wind stress. A 
siniilar result was obtained by Fissel et at. (1977). 
In this section we reexamine Hellennan's conclu- 
sion using the weather ship data and determine the 
effect of stability on the stress computation. 

The wind stress vector r can be determined from 
the bulk aerodynamic expressions 

7y ~ 

Here t, and Ty are the east -west components and 
north-south components of r, respectively, u, and 
Uy are the components of the wind; and Cu is the 
drag coefficient for momentum. The Liu et al. (1979) 
algorithm was used to compute Cu. The behavior 


Wind Stresses (dyn/cm^) for Jonuory 



Ffc. 8. Ctimaiologicsl wind stresses computed from the 
sampling method (abscissa) and a "direction-only" wind rose 
(ordinat'l. See also his. 3a. 
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Wind Stresses (dyn/ern^) for Jonuory 



Fig. 9. As in Fig. 8 except that the ordinate indicates the wind 
stress computed from the sampling method with neutral coef- 
ficients. 


of Co is similar to the behavior of Ch and Cu at wind 
speeds below 10 m s"' but increases significantly at 
wind speeds larger than 10 m s~'\ As described in 
Liu et al. (1979) momentum can be transferred by 
pressure forces acting on the air as it flows over the 
waves as well as by frictional stresses at the surface. 
The transfer by pressure forces tends to be strongly 
dependent on wind speed. However, molecular dif- 
fusion is the only mechanism available for the 
sensible and latent heat fluxes at the air-sea inter- 
face. Thus, the dependence of C« and Cg on wind 
speed tends to be less than the wind-speed depend- 
ence of Co for moderate to strong winds. 

Using the Liu et al. algorithm with the stability 
correction, we computed the climatological and 
Ty at each weather hip using the sampling method 
with a wind rose witn an average wind speed for each 
of the 16 direction categories. Fig. 8 shows that the 
"direction-only" wind rose consistently underesti- 
mates the wind stress when |tj] or |t„| is > 0.5 dyn 
cm"*. This relative error can be as large as 25% (see 
Tj. for Ship B). Thus, Hellerman's wind rose con- 
clusions for neutral coefficient wind stress also apply 
for stability corrected coefficients. Therefore, the 
wind stress should be computed either by the sam- 
pling method or by a wind rose with several speed 
categories for each direction. 

To clarify the imporance of stability we recom- 
puted the wind stress with the sampling method for 
a neutral Co. The results are shown for the clima- 
tological Tr and Ty in Fig. 9. The impoitance of 
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T, for January at Ship 0 



Fig. 10. Monthly averaged t , computed by the sampling 
method for stability corrected and neutral momentum transfer 
coefficients for each January at Ship D. 


Stability is small and affects the results by <109c 
except for small values of Tj. and r,. If we reexamine 
Fig. 6, we see that the effect of stability generally 
produces a smaller relative error for wind stress 
than for sensible and latent heat fluxes. 

The effect of stability on r^. and also was exam- 
ined on an interannual basis for each ship for January 
and July, and was found to be small in each case. 
This comparison is shown in Fig. 10 for the January 
T, at Ship D. It is apparent that for many applications 
the atmospheric surface boundary layer may be con- 
sidered to be neutral when computing wind stresses 
on time scales of one month or longer over the 
open ocean. 

The small errors caused by the use of neutral 
coefficients can be further reduced by substituting 
monthly mean values of (T, - T.) and (<j, - <?») 
into the transfer coefficient algorithms in place of the 
air-sea tempe' ature and humidity differences at each 
observation time. Errors of less than a few percent 
are typical of the results found at every ship position. 

8. Summary and discussion 

We have shown that monthly averaged values of 
middle-latitude and subtropical surface marine 
meteorological data can be used in bulk aerodynamic 
formulas to obtain the monthly averaged surface 
fluxes of sensible and latent heat. These results con- 
firm and extend the work of Budyko (1963) and his 
co-workers by including recent developments in 
boundary-layer similarity theory to give a more 


realistic representation of the physical processes 
which are important to air-sea transfers of sensible 
and latent heat. 

Our results also confirm that neutral transfer 
coefficients give reliable wind stress estimates over 
the ocean. If greater reliability is desired, this is 
easily achieved through the use of monthly aver- 
aged air-sea temperature and humidity differences 
for the determination of the drag coefficient. 

Similar results for heat flux and wind stress esti- 
mates should be obtained in tropical regions. Al- 
though this study does not include tropical marine 
data, our conclusions are valid for the small Bowen 
ratios (HILE) and the moderate to low wind speeds 
which characterize large portions of the tropical 
oceans (see. e.g.. Fig. 3b). 

We believe that the classical method for using 
wind roses and monthly averaged state variables 
to estimate heat fluxes and wind stresses has advan- 
tages over computing weighted sample means, and 
that these may prove valuable in describing and 
monitoring interannual and decadal climate fluc- 
tuations: 

1) Data processing requirements: The computa- 
tional effort required to process monthly-averaged 
data from ship-of-opportunity and ocean weather 
ships is much less than the effort needed to process 
all the individual surface repons. Therefore, the 
revision of heat flux estimates with an improved 
transfer coefficient or updated mean field is a rela- 
tively simple task with the classical method. 

2) Utilization of heterogeneous data: The classi- 
cal method may provide a simple way of combin- 
ing non-simultaneous data from a variety of observ- 
ing platforms. For example, satellite data does not 
yet provide accurate estimates of the air-sea tem- 
perature and humidity differences, but may soon 
give useful information on the surface wind speed 
over the oceans. To compute monthly averaged heat 
fluxes using the sampling method, it would be neces- 
sary to interpolate these data sets in space and time 
while at the same time combining the satellite data 
with conventional data sets. On the other hand, the 
use of the classical method relaxes the restriction 
that the observations of wind-speed, temperature 
and humidity fields be simultaneous and thus pro- 
vides a straightforward way to incorporate the 
satellite observations into the heat flux estimates. 

Some additional lines of investigation are sug- 
gested by the results of this study. First, the excel- 
lent agreement between the classical method and 
the sampling method for monthly averages raises 
the question of what is the shortest time interval 
for v'hich the conclusions hold. The study of Fissei 

at. (1977) begins to give an answer for a limited 
region over the North Pacific but this type of work 
should be clarified and should be extended to other 
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regions. Second, it is possible that the classical 
method may provide a more reliable estimate of the 
heat flux when the number of observations is small. 
Since both the covariance and the correlations of 
V, with either (T, - T,) or (q, - q,) are known to 
be relatively small, the nonlinearities in the flux esti- 
mates using the sampling method may lead to more 
statistiLsl uncertainly in estimating the long-term 
average heat fluxes than the classical estimators. 
Finally, the effects of area averaging on the con- 
clusions should be explored, particularly in regions 
of large gradients of surface parameters. 
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An atlas of the surface heat budget of the global ice-free ocean is 
presented. The atlas includes maps of the monthly mean distributions of 
the calculated heat budget components and of the basic variables used in 
the calculation. All the fields have a uniform spatial resolution and are 
presented on a regular 4 x5 grid. 

To calculate the heat budget componehts , bulk formulas parameteruing 
the heat flux mechanisms »ere applied to mohthly mean values of basic sur- 
face variables. The radiative fluxes were calculated using the bulk for- 
mulas documented by Budyko (1974), while the turbulent fluxes were 
calculated by the bulk aerodynamic method suggested by L1u et ul. (1979). 
Data sources for the basic variables were a global climatology of surface 
marine meteorological observations prepared by the National Climatic 
Center and the Berliand and Strokina (1930) cloudiness atlas. 

The results should be useful in ocean modeling and in a variety of 
diagnostic studies of the surface climate; all data presented are availabl 
in a computer compatible format. 
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1. Introduction 

The annual cycle of the earth's cll-ate clearly shows the Important 
role of the oceans. 3y storin, heat in the sumner season and releasinq i 
in the winter, and by transportinq heat from areas of heat gain in the t - 
pics to areas of heat loss at higher latitudes, the oceans moderate the 

Climatic state of the atmosphere. 

The energy sources for the ocean surface transport and storage are 
primarily the diabetic sources and sinks of heat at the ocean surface. 

These are the incoming solar radiation, the net loss by longwave radi tion 
and sensible and latent heat fluxes. Budyko was first to conduct a study 
of the spatial and temporal structure of the oceanic heat balance on glo- 
bal scale. Based upon bulk formulas obtained for this particular s u y. 
Budyko (1963) and his co-workers calculated the diabatic heating terms fro 
available observations of basic surface variables such as temperature, 

wind, humidity and sea level pressure. 

Since then many regional studies have been carried out by inves - 

tors using similar methods. Bunker and his co-workers (Bunker. 19 . 

Bunker and Worthington, 1976; and Bunker and Goldsmith. 1979) care u y 
constructed a relatively high resolution study of the ^ 

North and South Atlantic Oceans. Their study supported “ 

revealed by Budyko's study and offered better estimates on the str c 
of th heat exchange patterns especially over the Gulf Stream an 
Lth Atlantic. Bunker and his colleagues were also able to present local 

time series of the heat budget components. 

Hastenrath and Lamb (1977. 1978 and 1979) have used long te m x 

•latitude-longitude averages of surface marine ^ 

monthly high resolution climatology of the heat budget over part of t e 
Tropilal ocean. Hastenrath cabined these studies with others to Present 
the global aspects of the heat budget of the tropical ocean and atmospher , 

‘“*'‘r”rt’ly!^Welre at al. (1981) Studied the heat budget climatology over 
the tropical Pacific Ocean. They also addressed the problem “”'>'’"9 
errors inherent in using space and time averages of surface variables 

the bulk formulas (Weare and Strub, 1981). 
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Following a completely different approach Oort and Vender Haar (1976) 
addressed the question of the role of the oceans in the global heat 
budget and calculated meridional profiles of oceanic storage and northward 
heat transport over the northern hemisphere. 

While a combination of these studies provides a representation of 
oceanic heat budget, the results are not readily available for use in model 
studies of the general oceanic and atmospheric circulations. The only 
data source particularly designed for this purpose is that by Schutz and 
Gates (1971, 1972, 1973 and 1974) compiled from Budyko's (1963) study con- 
taining the months of January, April, July and October. The major goal of 
this study to obtain a similar presentation for all twelve calendar months 
for use by the modeling community. To achieve this purpose a readily 
available global climatology of surface marine observations with a resolu- 
tion of approximately 5x5 degrees was chosen and objectively checked and 
corrected. Bulk formulas were selected from the literature and a monthly 
climatology of the heat budget components was calculated over the entire 
ice free global ocean. These data and the climatology of basic surface 
variables are presented on a regular grid. 

2. Climatology of Basic Variables 

Budyko and his co-workers (Budyko, 1963) first utilized the so-called 
"classical method" for a comprehensive study of the global surface heat 
balance. Using monthly averaged atmospheric and oceanic data and bulk for- 
mulas they constructed the fields representing the spatial and temooral 
structure of the. various components of the balance. Their approach has 

been well documented by Budyko (1974). 

In this analysis we also use the classical method which allows the use 
of available surface climatological data. Two data sources were chosen, a 
global ocean climatology of marine observations prepared by the National 
Climatic Center (NCC), Asheville, North Carolina and an atlas of global 
cloudiness prepared by Berliand and Strokina (1980). In this section the 
data processing methods used in this study will be presented, followed by 
a brief discussion of the features revealed by the resulting fields. 
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2.1.1 The NCC climatology of marine surface ooservations 


The main source of data for this study was prepared by NCC based upon 
various sources of marine surface observations contained in the TDF-11, 
Surface Marine Observation, data collection!. These data have long served 
as a basis for the preparation of numerous marine surface atlases. They 
were also used in the regional ocean-surface heat budget studies referred 
to in the introduction. 

The specific climatology used here was prepared recently by NCC. The 
TDF-11 data were arranged into 1x1 degree Marsdin sub-squares and time 
averaged to obtain monthly mean values of the observed variables. These 
were then averaged into 5x5 degree sub-squares which were suitable for the 
purposes of this study. Table 1 describes the basic variables included in 
the 5x5 degree NCC climatology. 

The NCC 5x5 degree averages present climatological values which are 
assigned to the mean latitude and longitude of all observations within that 
square. This presentation may be thought of as an irregular grid whi-h 
changes from month-to-month. The transformation of the original climato- 
logy onto a regular grid was the first stage in this study. The gr'id cho- 
sen was that of the OSD Oceanic and Atmospheric General Circulation Models 
(hereafter referred .to as the OSU OGCM and A6CH). This is a rectangular 
latitude-longitude mesh with grid points located every 4 degrees in lati- 
tude and every 5 degrees in longitude^. Before starting the grid transfor- 
mation procedure, land, inland sea basins and the monthly sea-ice covered 
ocean points were marked so they could be masked out during the 

transf ormation^. 


iThe sources of data used in this collection, its contents and methods of 
quality control used during its assembly, are described in theTape Deck 
TDF-11 Manual published by the National Climatic Center, Asheville, North 
Carolina, 1968. 


2The latitude count starts at the south pole going northward and the 
longitude count at 180 longitude going eastward, there are 72 

points on each of the 46 latitude civ'cles, a total of 3312 grid points. 


3we based the land-ocean classification upon that of the OSU AGCM 
(Schlesinger and Gates, 1979). Sea ice boundary for each month was based 
upon data from Alexander and Mohlev (1976). 
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Table 1. Contents of NCC S'xS’’ Climatoloqy 


Location Information: (1) Marsden square number 

(2) S'xS* sub-square number 

(3) Mean latitude and lonqitude of 

observations 

Period; Month of the year 

Variable statistics: (1) Superstructure icing potential 

(2) Wind rose statistics - Monthly 
force and percent frequency of wind 
from eight points of compass and 
percent frequency of calms 

(3) Scalar variable statistics - monthly 
means, standard deviations and 
number of observations for the 
following variables: 

(a) wind speed 

(b) sea level pressure 

(c) air temperature 

(d) wet bulb tempe'-ature 

(e) dew point tenr -: ature 

(f) sea surface temperature 

(g) air - sea temperature difference 
(hj wave height 
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The qrld tr.insf ormatlon was achieved by usinq a Gaussian space 
averaging technique similar to what is known as the Barnes (1964) objective 
analysis scheme. The new grid point value is a weighted average of all NCC 
"grid" values falling within a radius of ^1100 km around It. The weights 
depend only on the distances between the new grid point and the mean loca- 
tion in each 5x5 degree Marsden sub-square entering the calculation. The 
Gaussian averaging technique is nearly mean conserving and acts like a spa- 
tial filter. It was designed to uniformly suppress features with wave- 
lengths shorter than '■1600 km. 

During the interpolation procedure, values based upon less than 25 
observations were immediately rejected, as were values located over the 
land and monthly sea-ice areas. A separate check was carried out on all 
values based upon 25 to 100 observations. If these deviated by more than a 
fixe-' amount from a value produced by a preliminary grid transformation, 
thej were marked ^or rejection and a final transformation was then per- 
formed without them^. 

For each variable field a corresponding data density field was also 
constructed. This was achieved by treating the number of observations of 
that variable as a field value and transforming it to the 4x5 degree grid. 
Then, the base 10 logarithm of the interpolated value was calculated at 
every grid point. These data density fields are used as data reliability 
indicators. 

The procedure described above was applied to 12 monthly sets of the 
following basic variables: sea level pressure, surface wind speed, sea 

surface temperature, sea - air temoerature difference, air specific humi- 
dity and sea - air humidity difference^. Air temperature fields were con- 
structed in a later stage from the monthly sea-surface temperature and sea 
- air temperature difference fields. 

4The maximum allowed differences between NCC values and the new grid values 
were 2.5 m sec~l for wind speed, 5 mb for pressure, S’C for SST, I’C for 
sea-air temperature difference, 3 g kg“^ for the humidity difference. NCC 
values were checked against the grid values after linearly interpolating 
from the grid to the NCC mean location of observations. 

^Air specific humidity and sea-air humidity difference were calculated 
prior to interpolating from the NCC dew point temperature and sea surface 
temoerature. Calculation was performed by usinq the Smithsonian (List, 
1953) formulas for saturation vapor pressure. 
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Examination of the resulting filtered data fields revealed relatively 
small wavelength features located in areas of low data density. These 
probably resulted from sampling errors due to the original area and time 
averaging procedures. In an effort to obtain a more reasonable represen- 
tation of the field we have applied additional smoothing to the fields in 
low data density areas. A special two dimensional Gaussian filter was 
designed which provided strongest smoothing in the zonal direction in lower 
data density areas. No smoothing was applied to areas where the data den- 
sity indicators were larger then 2.7 (500 observations). This procedure 
improved the representation esoecially in the southern hemisphere high 
latitudes where most of the problems occurred. It did result though in 
very small amplitude 2Ax features on the borderlines between areas where 
smoothing was applied and those where it was not applied. It is recom- 
mended that the reader bear in mind the daca density patterns when viewing 
or using the data fields. 

2.1.2 Cloudiness data 


The NCC climatology includes no cloudiness information and thus we had 
to use some other source of updated global cloudiness data. After com- 
paring two such sources. Miller and Feddes (1971) and Berliand and Strokina 
(1980), we decided to use the latter because it was free of some local 
deficiencies that show up in Miller and Feddfs's satellite study. Rerliand 
and Strokina* s data are also a more suitable inpiit for the bulk radiation 
formulas which were designed for use with sur.ace observations. 

The Berliand and Strokina data are based an a variety of sources, from 
regular ground observation networks over land and ocean to satellii-e obser- 
vatinns. These were all iniegrated and presented on a global rectangular 
5x‘^ oegree grid for eac^ cf the 12 months and are given both as figures and 
as tabulated values. The authors claim a better representation of the 
distributrion of cloudiness than achieved before, in particular over the 
ITCZ. 

For the purpose of this study the Berliand and Strokina data were 
transformed to a 4x5 degree grid using the basic Gaussian averaging tech- 
nique described in section 2.1.1. 
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2.1.3 Annual mean fields 


After completing the grid transformation and smoothing stages, an 
annual mean field was calculated for each of the variables. This was done 
by simply averaging 12 values at each of the annually ice-free ocean grid 
points. A data reliability indicator for the annual mean fields was 
constructed from its minimum monthly value at each grid point with the 
exception of the cloudiness fields which were not assigned data density 
indicators. 

2.2 DcL ^ L Z> 

Before discussing the main features in the resulting basic variable 
fields it is important to address the subject of the data quality. This is 
affected by errors of two categories: observation errors and sampling 

errors . 

The first category includes errors resulting from human factors and 
those resulting from instrumental limitations. Human errors most coimion in 
marine observations are those due to inadequate training and observational 
practice on behalf of the observers, and those introduced by the error 
susceptible communication links between the observer and the local climatic 
centers. Instrumental errors are commonly due to poor location of the 
instrument shelter or the sensors, and from the cruder instruments used on 
board merchant ships. These problems have been widely discussed in the 
open literature (see for example. Roll, 1965). The second category is 
associateo with non-uniform spatial and temporal distributions of the 
observations as well as from an insufficient number of observations in the 
statistical sense (resulting in very broad confidence intervals). These 
errors and their effrcts on determination of a long term climatology of 
basic and derived variables were addressed and examined by Weare and Strub 
(1981). As indicated by their study, these result in both random and bias 
errors. The bias errors are larger in areas of strong spatial gradients 
and where tne amplitudes of the diurnal and annual cycles are relatively 
large. The magnitudes of the biases also depend strongly on the number of 
observations used in the calculation of the climatological value. Weare 
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and Strub estimate that as far as basic variable climatology is concerned 
these errors are on the order of a few percent (and usually less than 10«). 

The automated data handling procedures used in this study result in 
a plausible representation of the basic variable fields. Large scale 
oceanic and atmospheric features are generally well resolved over most of 
the global ocean area. 

As indicated in section 2.1.1, field values over data-sparse areas 
should be viewed with some caution. In those areas some patterns show a 
considerable lack of time continuity. It should be noted, though, that only 
ocean areas south of about 50 $ are seriously plagued by these features. 

We believe however, that these data are still adequate for heat budget stu- 
dies on a global scale, in light of the uncertainties involved in esti- 
mating the heat fluxes from the semi-empirical bulk formulas. 

The data fields were compared with other available climatologies^ most 
of which based their marine climatology on the same TOF-11 data, though 
utilizing different data handling procedures and spatial resolutions. We 
found generally good agreement between these climatologies and ours. As 
a result of these comparisons we found that our automated data handling 
procedures performed well near the ocean - land boundaries. 


0 ^3 

£. * 0 


D^3aliSs^on o; 


bas'ic VGviahte fields tn z'ne atlas 


The basic variable fields constructed for the purpose of this heat 
budget study are presented in section 1 of the following atlas. Since the 
global distribution of these variables determines the global distribution 
of the heat budget components, it is important to examine the resulting 
fields for their agreement with other climatologies and for internal con- 
sistancy. This examination gives an indication of the inherent limitations 
in the quality and resolution of the heat balance components to be 
discussed later. 


^The climatologies used were; Schutz and Gates (1971, 

1974) Alexander and Mobley (1974), Hastenrath and l^mb 1977), and the 
Bunker and Goldsmith (1979) data. Some common marine atlases '^ere also 
used'^for comparisons in the various stages of the construction of the data 

handling routines. 
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2«3.1 Surface wind speed (Charts I.l to 1.13) 

The monthly-averaged surface wind speed fields generally agree with 
the monthly averaged sea level pressure (see Charts 1.92 to 1.104). Regions 
of weak winds are generally associated with the centers of the sub-tropical 
highs and regions of strong winds with the zones of steep pressure gra- 
dients near the mean positions of the middle latitude quasi-stationary low 
pressure cells. The agreement deteriorates in the Southern Hemisphere 
south of 50“S. The wind fields also reveal the trade wind belts as zones 
of local maxima on both sides of the ITC2 trough. The annual cycle asso- 
ciated with the Indian monsoon over the Arabian sea is also shown. 

2.3.2 Sea surface temperature (Charts 1,14 to 1.26) 

In these fields the patterns are usually more smooth and consistent 
probably due to a smaller variability and a larger number of observations. 
The warm and cold ocean currents are indicated by zones of sharp tem- 
perature gradients and deep narrow troughs and ridges penetrating poleward 
along the continental boundaries. Strongest temperature gradients are 
associated with the Gulf and Kuroshio currents. These gradients show a 
Hiarked annual cycle associated with the reduced horizontal temperature 
contrast during the summer months. Another belt of sharp temperature gra- 
dient is located over the Antarctic Ocean basin surrounding the Antarctic 
continent. It is more difficult to define seasonal variations in this 
region because of the lower data densities involved. Large amplitude 
equatorial features are seen in the eastern oceans. Here the annual cycle 
is also well represented by the shift of the equatorial maximum. The whole 
equatorial ridge especially the zone of warmest water over the Indian and 
Central West Pacific Oceans oscillates in position throughout the. year. 
Finally, smaller scale features like the southward penetration of the 
Labrador current or the northward penetration of the Alaska current are 
also well defined. 
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2.3.3 Sea - air temperature difference (Charts 1.27 to 1.39) 

The outstandinq patterns in these fields are those associated with the 
large contrast between the winter continental air masses and the warm ocean 
surfaces on the western sides of the northern hemisphere ocean basins. 

Here the instability of the atmospheric boundary layer is maximum in 
December and is accompanied by strong horizontal temperature difference 
gradients. These figures reveal the stabilizing effects of the cold water 
currents to the north of the Gulf and Kuroshio. This is especially strong 
during summer when a stable air mass is created over the Labrador and 
Kamchatka currents. In general, the atmospheric boundary layer over the 
suirmer hemisphere is more stable than that of the winter hemisphere. 

2.3.4 Air temperature (Charts 1.40 to 1.52) 

These fields were constructed from the sea surface temperature and 
sea - air temperature difference fields and as such contain no additional 
information. 

2.3.5 Air specific humidity (Charts 1.53 to 1.65) 

The specific humidity of the air is strongly dependent on its tem- 
perature. This is shown by the strong similarity between these fields and 
the air and sea surface temperature fields. The spatial gradient in the 
humidity fields though smaller than in the temperature fields, are still 
directed from polar to equatorial regions. The strong zonal structure of 
these fields is, slightly modified by troughs and ridges in the north-south 
direction associated with the oceanic circulations. 

2.3.6 Sea surface - air specific humidity difference (Charts 1.&6 to 1.73), 

This variable is sometimes referred to as the evaporation potential. 

It is generally largest in the subtropical zones where subsiding dry air 
masses overlay relatively warm water surfaces. These zones usually extend 
over the southern branches of the warm Gulf and Kuroshio currents, except 
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during the northern hemispheric suimier months when the surface boundary 
layer is relatively stable and its vertical humidity gradient less steep. 
The ITCZ minimum is enhanced over the cold upwelling regions of the eastern 
Atlantic and Pacific Oceans. These minima become especially sharp during 
the northern hemisphere surmier months. Also noticeable are the strong 
horizontal gradients in the evaporation potential on the poleward boun- 
daries of the middle latitude warm ocean currents. Some local features of 
interest are the maximum over the western coast of Central America per- 
sistent throughout the year which is associated with the local maximum in 
sea surface temperatures and the sea-sonal changes over the Arabian sea 
associated with the monsoon circulation. 

2.3.7 Cloudiness (Charts 1.79 to 1.91) 

These fields generally show maxima over the high latitudes and minima 
in the subtropics. The ITCZ is indicated by a local shallow ridge over the 
tropics, except for the tropical western Pacific where the cloudiness 
broadens to cover a large ocean area. T.ones of marine fog and stratocumu- 
lus prevailing over the cold ■^astern ocean basins are well depicted. These 
also show stronger seasonal dependence than most other features. Overall, 
the cloudiness fields have a smoother appearance than ' c basic variables 
from the NCC climatology. 

2.3.8 Sea level pressure (Charts 1.92 to I.IQA) 

The monthly averaged sea level pressure fields present the familiar 
patterns of the oceanic subtropical highs and the seasonal changes in their 
strength and location. In the Northern Hemisphere the Aleutian and 
Icelandic lows and their seasonal changes are also well depicted. Also 
indicated are the seasonal changes in sea level pressure over the Arabian 
Sea wh.ch are associated with the Indian Monsoon. In these coarse resolu- 
tion fields the ITCZ appears as a broad and shallow pressure trough located 
over the tropics and migrating northward or southward with the seasons. 

Not so well represented are the pressure patterns over the southern 
hemisphere high latitudes. This is indicated primarily by lack of time 


continuity of some of the patterns. These patterns are associate! with low 
data density and are usually confined to south of 50 S. 

3. Ocean Heat Budget Climatology 

The local monthly-averaged energy balance of the ocean's upper layer 
can be expressed as follows: 

S = -OivF - h’w + Fs(l - - Fi - Fl£ - Fh 

where S Is the layer-averaged rate of change of heat storage per unit area 
of the layer (storage term), OlvF is the layer-averaged divergence of heat 
transported by the surface horizontal currents, and is the downward heat 
transport through the lower boundary of the layer. The remaining terms are 
the diabatic sources and sinks of heat at the ocean surface: Fj is the 

incoming solar flux at the ocean surface reduced by the monthly mean sur- 
face albedo Fj is the net upward long-wave radiative flux, F(_£ the 
upward latent heat flux, and Fh the upward sensible heat flux at the ocean- 
atmosphere interface. In constructing the monthly estimates of the diaba- 
tic terms, the classical approach first used by Budyko (1963) was taken, 
i.e., the monthly averaged fluxes were calculated from the monthly averaged 
values of the basic variables. The calculations were carried out on the 
4x5 degree grid using the previously described climatology. The heat 
fluxes were calculated separ .tely at each grid point. All the fields that 
were based on variables taken from the NCC climatology were also assigned a 
data quality indicator. This is defined as the local minimum of all indi- 
cators associated with the input variables. 

Estimating the diabatic heating sources will enable us to calculate 
the sum of the monthly storage and the flux divergence terms in the surface 
layer as a residual. When averaged annually and globally this sum should 
vanish, and an approximate balance should be achieved between the diabatic 
fluxes if we assume that the contribution of Fyy is negligible. 


ORIG’N^'- 
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3.1 Estimating the components of the heat budget 
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3.1.1 Incoming solar radiation 

There are several methods presented in the literature which are 
intended to provide an estimate of the solar flux incident upon a horizon- 
tal surface at sea level. AH of these consist of two parts: an estimate 

of the monthly averaged maximum available direct plus diffuse solar flux at 
sea level, which is dependent on latitude and season, and a factor which 
includes the effects of cloudiness. Of all the bulk formulas used for cli- 
matological studies, only that suggested by Berliand (1960) was designed 
for a comprehensive global study (Budyko, 1963), under the constraint of 
achieving a global heat balance. 

The Berliand algorithm is described rigorously by Kondratyev (1969) 
and Budyko (1974). It is easy to apply and combines both parts of the 
estimate in a consistent manner. Although its clear sky formulation and 
cloudiness factor performances have been criticized by Reed (1977) and 
others based on observations at specific times and locations, we believe 
one should judge its performance as a whole over a wide range of conditions 
as indicated by Budyko (1974). In their tropical Atlantic study Hastenrath 
and Lamb (1978) used a different method to estimate Fj but their estimates 
were no more than a few percent away from those obtained by Bunker (1976) 
and Bunker and Goldsmith (1979) using Berliand's algorithm. Reed's (1977) 
suggested method was used by Weare and ‘Strub (1981) with results 10% larger 

than those calculated here. 

Berliand's formula as presented by Budyko (1974) is 

Fs = Qotl - (a bn)n] 

where Qo is the monthly averaged direct and diffuse maximum solar flux and 
is a function of latitude, n is the monthly averaged cloudiness (in 
fractions), b is a constant and a is a coefficient varying with latitude. The 
monthly values of Qo were compiled by Berliand from observations from 
various latitudinal zones and are given by Budyko (1974, Table 3). The 
coefficients a and b were also calculated by Berliand and are given in 
Budyko's book (a values are presented in Table 4 and b = 0.38). 
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To assess the amount of solar heating of the ocean surface one has to 
subtract the reflected part of the incoming solar flux. For that purpose 
Payne's (1972) albedo table for the Atlantic was used globally as suggested 
by Kukla and Robinson (1980). 

3.1.2 Net long-wave radiation 


Here again there is a large diversity and controvery in the literature 
as to hew to construct a reliable estimate of the local flux. Most of the 
methods agree closely with that originally suggested by Brunt (1932) and 
later modified somewhat by Berliand and Berliand (1952). This method, as 
presented in Budyko (1974), was the only one used for a global comprehensive 
heat budget study. Again the formula consists of two main parts: an 

estimate for the cloud-free net IR lost per unit area and time by the ocean 
surface, and a cloudiness factor. The formulation of the cloudiness effect 
is the most controversial aspect and seems to be the main source for the 
differences between estimates obtained by different investigators. 

Berliand and Berliand's formulation is a linearized version of the 
difference between the grey body fluxes of the ocean (upward) and 
atmosphere (downward). The downward atmospheric flux should be affected 
mainly by the amount of water vapor in the lower layers and the amount and 
type of cloudiness. In Berliand's formulation the net flux is given by: 


Fi = £cjTa(p - q*''e)(l - cn"^) + 4eaTa(Ts - T^) 


where T^ and T 3 are the monthly mean air and sea surface temperatures, e is 
the monthly averaged water vapor pressure at anemometer level (in mb), and 
n is the mean cloudiness (in fractions). £ is the emissivity of water 
(=0.97), CT is the Stefan-Boltzmann constant, p anc q are constants 
(according to Berliand, ^o:ual to 0.39 and 0.05, respectively), c is a 
coefficient depending on latitude in order to account for the varying 
radiative properties cf clouds at different latitudes; Berliand assumed a 
linear dependence on latitude? ind his values were used here (Budyko, 1974, 
Table 9). t is a constant and was set equal to 2 in this study. 
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Simpson and Paulson (1979) found < 50 od agreement between estimates 
based on Berliand's form- of Brunt's formula and observations under clear 
sky conditions. Weare and Strub (1981) used a different cloudiness factor 
while basically not changing other parts of Berliand's expression. 
Comparing their final results with ours we found that their annually 
averaged Fj values are smaller by 20 to 40?i. Here again we decided to 
apply Berliand's formulation as a whole, for the same reasons as before, 
namely, comprehensiveness and design for global and long term studies. 

3.1.3 Sensible and latent heat fluxes 

The latent and sensible heat flux estimates were based upon the bulk 
aerodynamic formulas 


Fh = i3CpCHVa(Ts - Ta) 


and 


Fle = CEVaCqs - Qa) 

where o is the density of air, Cp the isobaric specific heat of air, and L 
the latent heat of evaporation^. V^, Tj and qa are the monthly averaged 
wind speed, temperature and specific humidity of air at anemometer level. 

Tj is the sea surface temperature and Qj the humidity of the layer of air 
close to the sea surface, assumed saturated at temperature T 5 . Ch and 
Ce are the transfer coefficients for sensible and latent heut; these were 
calculated using the iterative method suggested by Liu et al. (1979), using 
the observed T 5 for the "skin" temperature of the ocean. This method was 
tested by Esbensen and Reynolds (1981). 

^Here 0 was calculated locally from the monthly mean pressure and virtual 
temperature. Cp and L were taken as constants equal to their value at 0 C. 
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3.1.4 Some general remarks 


The bulk formulas used to calculate the diabatic fluxes are parame- 
terizations of the physical mechanisms involved in the heat transfer pro- 
cess. As such they are approximations to the true flux values. The 
approximation errors are believed to be in the range of 10 to 20?». Such 
errors are larger than bias errors arising from irregular sampling of basic 
variables in time or space, as was shown by Weare and Strub (1981). Their 
studies indicate that the magnitudes of the bias errors are strongly depen- 
dent on the number of observations of the basic variables and the geo- 
graphic location (through the local spatial gradients of the variables and 
the magnitudes of their diurnal and annual cycles). Weare and Strub esti- 
mate that such sampling errors as reflected in the heat flux estimates are 
about lOlJ. 

3.2 Discussion of I'esulzs 

The result of the heat flux calculations are presented in part 2 of 
the atlas. These include monthly sets of the spatial distributions of 
the fluxes and the calculated residual identified as the net downward heat 
flux. For each month a net downward radiative flux (Fs(l - ^) - Fj) is 
also presented. Based upon these monthly fields a set of zonal mean vs. 
time sections of the spatial distribution of the local annual mean (only 
for annually ice free grid points) are presented for each of the com- 
ponents. Finally, the annually averaged meridional profiles of the com- 
ponents are presented based on the monthly zonal mean values for each of 
the components. 

3.2.1 Available solar heat flux (Charts 2.1 to 2.13) 

These fields present the distribution of the net solar flux (including 
the albedo effect) at the ocean surface. This is dependent upon the lati- 
tudinal maximum solar zenith angle and the distribution of cloudiness, the 
latter determining its longitudinal structure. Thus we find maxima in 
available solar flux over the relatively cloud-free subtropics with the 
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annual maximum occuring during the northern hemisphere summer over the 
western coast of North Africa. The ITCZ can be identified as a narrow zone 
of local minimum but is sometimes obscurred by larger scale features. The 
seasonal changes result mainly from changes in the solar zenith angle and 
are shown as an increase in the flux values over the middle latitudes 
during the sumner season and a decrease during winter. Seasonal changes 
also occur in the relative strength of the subtropical maximum in both 
sides of the equator. 

3-2.2 Net upward long-wave flux (Charts 2.14 to 2.261 

Over the warm and relatively cloud-free subtropics the net IR fluxes 
reach maximum values. Other local maxima occur over the warm ocean 
currents in winter, where the ocean loses more than it gains in longwave 
radiation. The ITCZ is a zone of local minimum since the relatively cloudy 
sky and colder ocean surfaces result in increased downward radiation and 
reduced upward radiation simultaneously. The strongest annual cycle is 
observed over the warm ocean currents. Here factors like the decrease in 
ocean - air temperature contrast and a more humid atmospheric boundary 
layer result in a marked decrease in the net IR flux leaving the ocean sur- 
face during the sutimer months. These changes are strongest over the Gulf 
and Kuroshio currents but are also noticeable over the Antarctic Ocean. 

The seasonal changes over the tropical oceans are small. 

3-2.3 Net downward radiative flux (Charts 2.27 to 2.39) 

These are presented as another means of observing the contribution of 
radiative fluxes to the heat budget. 

3-2.4 Latent heat flux (Charts 2.40 to 2,52) 


The flux of latent heat associated with the evaporation process pre- 
sents the largest heat sink for the ocean surface layer. It is a strong 
function of the evaporation potential Ps-Qa l^ut also depends on the wind 
speed and the stability of the lower part of the atmospheric boundary 
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layer. The major patterns and their annual cycles closely resemble those 
of the sea - air humidity difference fields (Charts 2.79 to 2.91). The 
major maxima occur over the subtropics and extend northward over the Gulf 
and Kuroshio currents. These ridges weaken during the northern hemisphere 
summer. Fluxes are relatively weak over most of the low latitude oceans 
because of the stability of the boundary layer and the weak surface winds. 
The major minima occur over the regions of upwelling in the eastern 
equatorial Pacific and Atlantic oceans. A local maximum exists over the 
Arabian sea during the summer Monsoon months. The largest annual changes 
occur over the warm ocean currents of the Northern Hemisphere. 

3.2.5 Sensible heat flux (Charts 2.53 to Z.55) 

The sensible heat flux is a strong function of the sea - air tem- 
perature difference and as such their spatial and temooral distributions 
are similar. In the global average sensible heat flux is the weakest heat 
sink for the surface layer. Large fluxes occur regionally over the Gulf 
and Kuroshio currents during the northern hemisphere winter and near the 
southern tip of Africa out over most ocean areas these fluxes are very weak 
and show little coherent spatial structure. 

3.2.6 Net downward heat flux (Charts 2.66 to 2.78) 

This residual of the diabatic heating terms is mostly influenced by 
the patterns of solar radiation and latent heat fluxes. Large heat losses 
are found over the warm water currents while strong gains exist over the 
cooler surfaces of the tropical oceans. On the average both hemispheres 
lose heat during the winter months and gain during surrmer. The annual mean 
field presents the part of the residual available for transport on an 
annual basis. This field of course shows the tropics as a heat source and 
the high latitudes as heat sinks. The largest sinks are found over the 
western boundary currents and the largest sources are found in the regions 
of :^rong equatorial and coastal upwelling. It is interesting to note that 
this study shows that on the average the southern at 1 antic ocean does not 
act like a sink but rather like a weak source of heat. 
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3.2.7 Zonal means vs. time sections (Charts 2.79 to 2.S4) 

These figures illustrate the seasonal changes of the heat budget com- 
ponents. Zonal averaging was done over ice-free points only. 

3.3 Conaluiiing 

To summarize the study we present the meridional profile of the ann- 
nually averaged net downward heat flux and its global mean (Fig. 1). These 
lead to an average residual heat storage rate of 5.37 Wm"2. This value 
indicates that the results comply with our preconceived notion of the 
existence of a global balance. Here we should remember that only ice-free 
points are included in this study^. It is most likely that the ocean heat 
loss to the ice covered ocean would tend to balance part of the computed 
storage rate over the open ocean. In any case the residual amount of heat 
distributed over a 1 km deep ocean would result in an annual temperature 
increase of -0.05 C. 

The meridional profiles of the heat budget components are presented in 
Fig. 1 and Table 2. The shape of the net downward flux profile agrees 
cuite well with that presented by Budyko (1974, Table 19) and Hastenrath 
(1980). Table 3 presents the total ocean area vs. the ice-free ocean area 
in this data set. 

Finally, we t.ote that the boundary layer fluxes of sensible and latent 
heat are systematically smaller than those calculated by Budyko (1963) and 
Bunker (1976). However, we believe that the transfer coefficients adopted 
for this study give a more realistic description of the sensible and latent 
heat fluxes with changes in wind speed and boundary layer stability. It is 
also encouraging to note that approximate annual heat balance can be shown 
over the global oceans with the new formulation. 

^The annual-average zonal mean values were calculated as follows: 

i: i: 

Z * m»l ‘ m^l 

where Z^ is the monthly zonal average taken over the ice-free ocean points 
only and Nn, is the monthly number of ice-free ocean points in the same 
latitude belt. 
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Meridional profiles of the components of the ice-free 
ocean annual heat budget (in W m"-). and Fq stand 
for net radiative flux and net downward heat flux and the 
overbar for an average. Other sytnbols have the same 
meaning as in text. 
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Table 2. 10® zoO'il means of the components of the ice-free ocean 

annual heat budget (in W m“2). Fq is the net downward heat 
flux. Other symbols have the same meaning as in text. 


^ ^ ^ Ps(l- *w) ^ ^ 


80N-90N 


WWW 

WWW 

WWW 

WWW 

WWW 

70N-80N 

29 

50 

46 

31 

35 

-44 

60N-70N 

23 

51 

51 

39 

38 

-36 

50N-60N 

16 

48 

40 

33 

35 

-29 

40N-50N 

12 

61 

51 

108 

57 

-16 

30N-40N 

14 

113 

58 

160 

102 

-25 

20N-30N 

9 

131 

57 

193 

136 

- 4 

10N-20N 

4 

131 

57 

195 

146 

12 

0 -ION 

4 

109 

46 

191 

145 

32 

lOS- 0 

2 

106 

43 

199 

151 

43 

20S-10S 

3 

125 

52 

194 

142 

14 

30S-20S 

6 

118 

58 

181 

123 

1 

40S-30S 

6 

96 

59 

157 

98 

- 4 

50S-40S 

2 

62 

53 

114 

61 

- 3 

60S-50S 

2 

43 

43 

79 

36 

- 9 

70S -60S 

5 

34 

39 

SO 

41 

2 

30S-70S 

6 

26 

38 

74 

36 

4 

90S-80S 

*** 

★ ★★ 

WWW 


WWW 

WWW 

GLOBAL 

6 

98 

51 

160 

109 

5 



-B22- 


Table 3: Total ocean area (1) and annually averaqed ice-free ocean 

area (2) in 10^ latitude belts in units of 


( 1 ) ( 2 ) 


80N-90N 

70N-80N 

60N-70N 

50N-60N 

40N-50N 

30N-40N 

20N-30N 

10N-20N 

0 -ION 
lOS- 0 
20S-10S 
30S-20S 
40S-30S 
50S-40S 
60S-50S 
70S-60S 
80S-70S 
90S-80S 

GLOBAL 


3.207 

7.074 

4.491 

10.214 

14.065 

13.796 

24.083 

30.743 

34.337 

34.338 
33.356 
30.634 
31.991 
30.735 
25.424 
17.236 

4.126 

0.000 

355.351 


0.000 

1.881 

3.512 

9.620 

13.955 

13.796 

24.083 

30.743 

34.837 

34.338 

33.356 

30.634 

31.991 

30.735 

23.982 

7.364 

0.102 

0.000 

329.931 
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4, Summary 

In this study the global ocean heat budget was estimated using bulk 
formulas and available climatologies of surface observations. The objec- 
tive approach to the analysis renders such a global study simple and easy 
to update if additional data or improved formulas become available. 

Recent developments in boundary layer theory were used for calculating 
the boundary fluxes of sensible and latent heat. This method leads to 
lower estimates of the fluxes when compared with previous studies, but an 

approximate global balance was achieved. 

The bulk radiative formulas used were those suggested by Budyko and 
his colleagues during the preparation of their global study (Budyko, 1963). 
These have been seriously criticized in the open literature but none of the 
suggested alternatives have been subjected to a test or achieving a global 
balance against reasonable values of the turbulent fluxes. We believe that 
the radiative fluxes in this atlas are reasonable and that the next advance 
in determining the surface heat balance will not come from the improvement 
of the bulk radiative formulas for surface observations. Rather it remains 
the task of future studies to offer a direct and more accurate estimate of 
these fluxes through the full utilization of satellite information. 

The fields of basic variables and heat budget components constructed 
during this study provide an easily accessible data source for scientists 
interested in boundary conditions for atmospheric and oceanic general cir- 
culations models or in evaluating model performance. 

The complete data set presented in this report is available on a com- 
puter compatible tape from the Climatic Research Institute, Oregon State 
university, Corvallis, Oregon, 97331. 
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MONTHLY VARIATION OF THE ZONALLY-AVERAGED HEAT BUDGET COMPONENTS 
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